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Abstract: We explore the possibility to observe the effects of electron neutrinos from 
past galactic supernovae, through a geochemical measurement of the amount of Tech- 
netium 97 produced by neutrino-induced reactions in a Molybdenum ore. The calcu- 
lations we present take into account the recent advances in our knowledge of neutrino 
interactions, of neutrino oscillations inside a supernova, of the solar neutrino flux at 
Earth and of possible failed supernovae. The predicted Technetium 97 abundance is of 
the order of 10 7 atoms per 10 kilotons of ore, which is close to the current geochemical 
experimental sensitivity. Of this, ~ 10 — 20% is from supernovae. Considering the 
comparable size of uncertainties, more precision in the modeling of neutrino fluxes as 
well as of neutrino cross sections is required for a meaningful measurement. 
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1. Introduction 

The question of the rate of core-collapse supernovae in our galaxy is as interesting - 
in connection with the history of star formation, the origin of elements and climate 
evolution - as it is difficult. 

Next to more traditional approaches - based on astronomical searches, observa- 
tion of radioactive decays and historical records - an interesting combination of neu- 
trino physics, nuclear physics and geochemistry offers a way to count past supernovae 
and learn about their neutrino output. The idea, first proposed in 1988 by Haxton 
and Johnson (HJ) [1], relies on the fact that energetic electron neutrinos can produce 
Technetium 97 ( 97 Tc) from natural Molybdenum ( 98 Mo and 97 Mo) via the reactions 
98 Mo(^ e , e - ri) 97 Tc and the 97 Mo(z/ e , e~) 97 Tc. The lifetime of 97 Tc is 2.6 ■ 10 6 years, 
which eliminates the possibility of the 97 Tc found on earth to be primordial. There- 
fore presence of 97 Tc in a Molybdenum ore is a result of either neutrino or cosmic-ray 
reactions. Any excess of 97 Tc, over the expected contribution of solar neutrinos and 
of cosmic-rays, is an historical imprint of nearby supernova explosions that happened 
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before any human technology was there to study them! This suggests the possibility 
for an "archaeological" quest of the neutrinos from past galactic supernovae, through 
a geochemical measurement. 

In spite of the challenges and uncertainties involved, Haxton and Johnson came to 
positive conclusions: the neutrino-induced production of 97 Tc in Molybdenum rocks, 
such as the Henderson mine in Colorado, was estimated to happen at an average of 
1.57 • 1CT 36 captures per second per target atom, about 40% of the background due to 
solar neutrinos, and corresponding to a detectable abundance of 97 Tc at the present 
time. It was understood that the background due to cosmic rays could, in principle, be 
reduced to acceptable levels by excavating deeply enough and/or measuring the 97 Tc 
attenuation rate as a function of the excavation depth. 

In 1988 the idea was ahead of its time, in a field mostly focused on the disappearing 
solar and atmospheric neutrinos. Therefore it remained restricted to a small community 
[2]. Initial experimental efforts, led by K. Wolfsberg (LANL), lacked the funds necessary 
to overcome technical difficulties, and were soon abandoned [3, 4, 5]. 

Today the field is more mature: after solving the solar and atmospheric anomalies 
and discovering neutrino oscillations, research has focused more on supernova neutrinos, 
whereas supernova models have become sufficiently predictive, a geochemical measure- 
ment of 97 Tc production could constrain the galactic supernova rate (the extragalactic 
contribution is negligible), for a given model of neutrino spectra and luminosities from 
an individual supernova, without the extinction effects typical of astronomy. Notably, 
even dark neutrino sources - e.g., those that give no luminous signal due to rapid for- 
mation of a black hole - would be counted. The capability to trace back to ancient 
supernovae in our galaxy would be an unique complement to the activity of real time, 
large scale liquid detectors. With the masses of 0.1-1 Mt envisioned for the near future 
[6, 7], these detectors should see the diffuse flux of neutrinos that continuously reaches 
the Earth from all the supernovae of the universe and that has a substantial cosmo- 
logical component (see e.g. [8]). Finally, the 97 Tc data could be the only one available 
on electron neutrinos from supernovae! They would add complementary information 
compared to the v e data from SN1987A [9, 10, 11]. It is fascinating that this infor- 
mation, so important to have a complete picture of core-collapse, is there, buried in 
Molybdenum rocks, only waiting to be deciphered. 

In the present paper we perform detailed calculations of the expected amount of 
Technetium 97 produced by the neutrino-induced reactions in Molybednum ore, due to 
both solar and supernova neutrinos. Our predictions are based on the best available 
results on the solar neutrino fluxes and on the recent developments in the study of 
neutrino oscillations inside the star. To estimate uncertainties, we also consider a 
range of values for the galactic supernova rate as well as for the neutrino fluxes at the 
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neutrino-sphere, in accordance with supernova modeling. The structure of the paper 
is as follows. In Section II and III we give generalities on the supernova galactic rate 
and on the solar and supernova neutrino fluxes. Section IV describes calculations of 
the neutrino-nucleus cross sections. The results are summarized in Section V, while 
conclusions are in Section VI. 

2. Super novae in our galaxy: rates and distances 

The rate of core collapse supernovae (SNe) in the Milky Way is known only up to its 
order of magnitude. The rate R sn ~ 0.09 yr" 1 used by HJ, is practically excluded by 
the lower values measured recently (table 1). Scaling from external galaxies and gamma 
rays from galactic 26 A1 independently favor a rate of R sn ~ 0.02 yr -1 as central value, 
and up to R sn ~ 0.05 yr _1 at 3a. These are consistent with the much more uncertain 
estimate from historical records of galactic SNe, which favors a higher rate. The upper 
limit from the non-observation of supernova neutrino bursts is rather loose, but still 
of interest for its being extinction-free and sensitive to possible failed supernovae, as 
mentioned in sec. 1. 



Table 1: Estimated rate of galactic core-collapse SNe per century, from [12], with permission 
o f the author. 



Method 


Rate 


Authors 


Refs. 


Scaling from external galaxies 


2.5 ±0.9 


van den Bergh 


[13, 14] 






& McClure (1994) 






1.8 ± 1.2 


Cappellaro & Turatto 


[15, 16] 






(2000) 




Gamma-rays from galactic 26 Al 


1.9 ±1.1 


Diehl et al. (2006) 


[14] 


Historical galactic SNe (all types) 


5.7 ± 1.7 


Strom (1994) 


[17] 




3.9 ± 1.7 


Tammann et al. (1994) 


[18] 


No neutrino burst in 25 years" 


< 9.2 


Alekseev & Alekseeva 


[19] 




(90% CL) 


(2002) 





"The limit of Ref. [19] is scaled to 25 years of neutrino sky coverage. 



Here we keep R sn as a free normalization; results for transition rates will be given 
for R sn ~ 0.01 yr _1 , as they can be easily rescaled for other values. The final results 
for the 97 Tc abundance (table 6) refer to the more realistic R sn ~ 0.03 yr _1 , for easier 
comparison with previous literature on supernova neutrinos where this rate is commonly 
used. 
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To describe the spatial distribution of supernovae in our galaxy we use cylindric 
coordinates (r, z, 9) centered at the galactic center, which is at d & ~ 8.5 kpc of distance 
from us. We assume a uniform distribution in 9 and adopt the following function for 
the distribution in r and z: 



n cc (r) oc r^exp(— r/u) x 



0.79 e ~ {z/212 pc)2 + 0.21 e ~ (2/636 pc)2 



(2.1) 



which fits observations of neutron stars with the parameters £ = 4 and u = 1.25 kpc 
[20, 21, 22], and is illustrated in more detail in ref. [23]. 




Figure 1: The distribution with the distance from Earth of a Monte Carlo-generated 
population of 10 5 supernovae. The distribution in Eq. (2.1) was used, with £ = 4 and 
u = 1.25 kpc. 

With this, we calculate the distribution of supernovae with the distance from Earth, 
using both numerical integration and the Monte Carlo Method (fig. 1). We find that the 
neutrino flux at Earth corresponds to the flux one would obtain if all the supernovae 
were at the effective distance d e ff = 4.26 kpc, which we use throughout the paper. 
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With the alternate set of parameters given in [23] (£ = 2.35, u = 1.528 kpc) we get a 
value only slightly different, d e // = 4.6 kpc. 

The Monte Carlo approach was used to investigate possible variations in the results 
due to one single supernova very close to the Earth, d ^ 100 pc or so. Our results agree 
with those in HJ on how such variations can be safely neglected: only very few (~ 2 
in 10 5 ) events satisfy this condition, therefore the uncertainty due to the fluctuation of 
supernova positions in the galaxy turns out to be small compared to other uncertainties 
of different nature (total neutrino luminosity per supernova, etc.). 

We have not investigated effects of local deviations from the spatial distribution 
(2.1). Nguyen and Johnson have studied this possibility [2], pointing out that a par- 
ticularly high supernova rate in the Sco-Cen region of the galaxy might result in an 
increase of the 98 Mo(f e , e~n) 97 Tc production rate, to a level equivalent to the whole 
galactic supernova neutrino production with R sn = 0.01 yr" 1 rate. The full investiga- 
tion of this possibility is beyond the scope of our paper; still one can effectively include 
it here as a further uncertainty in the rate R sn , which is a free parameter for us. 

3. Solar and supernova neutrino fluxes 
3.1 Neutrinos from supernovae 

A core collapse supernova is an extremely powerful neutrino source. It releases about 
3 • 10 53 ergs of energy within ~10 seconds in neutrinos and antineutrinos of all flavors 
before the star explodes, leaving behind a neutron star or, in rare cases, a black hole. 

The spectrum of the neutrinos of a given flavor w at production in the star is 
approximately thermal, and commonly described by the form [24]: 



with average energies Eq w in the range of 10 — 20 MeV and with the muon and tau 
species (collectively called v x from here on) being harder than the electron ones due 
to their weaker (neutral current only) coupling to matter. L w is the (time-integrated) 
luminosity in each flavor, L w ~ 5 • 10 52 ergs; a w is a parameter describing the shape of 
the spectrum, a w ~ 2 — 5 [24] , with larger a w corresponding to narrower spectrum. 

For the purpose of studying the dependence of our results on the original neutrino 
fluxes, we choose three sets of parameters (Table 2) that are inspired by current nu- 
merical calculations [25, 24, 26] (see also [8] for a summary of those) and that can be 
considered as representative of natural and extreme cases. The cases of an especially 




(3.1) 



5 



high or low v e flux at Earth are modeled mostly by varying the average energy and lu- 
minosity of the original v x flux. This is justified by the fact that this flux is responsible 
for at least ~ 70% of the v e flux in a detector, as will be clarified below. 





best 


natural 


worst 


(E e ,E x )/MeV 


(13,22) 


(12,18) 


(10,16) 


{L e , L x )/ L 


(1,2) 


(1,1) 


(1, 0.5) 


(oc e ,a x ) 


(3.5,2.5) 


(3.5,2.5) 


(3.5,2.5) 



Table 2: The three sets of input parameters for the neutrino fluxes used in this paper. Here 
L = 0.5 • 10 53 ergs. 

As they propagate from the production point to a detection point on Earth, the 
neutrinos undergo flavor conversion (oscillations). Therefore, up to a geometric factor 
I /And 2 , with d the distance star-Earth, the flux of u e in a detector at Earth depends 
on the fluxes at production as: 

F e = pF° + (l-p)F° x , (3.2) 

with p the u e survival probability. For a given density profile of the star p is determined 
by number of physical phenomena, namely: (i) neutrino -neutrino coherent scattering, 
which induces complicated non-linear effects [27, 28, 29, 30, 31, 32, 33, 34, 35, 36] (ii) 
two MSW resonances produced by neutrino-electron scattering [37, 38, 39] and (iii) 
oscillations inside the Earth [38, 40, 41]. As a result p is a function of the neutrino 
energy, the mixing angle 9^ and the mass hierarchy (ordering) of the neutrino mass 
spectrum, i.e., the sign of the mass squared splitting Am^ (see e.g., [38]). For the 
normal mass hierarchy, Am^ > 0, p varies in the interval 

p = - sin 2 # 12 ~ - 0.31 , (3.3) 

as sin 2 #13 increases in the interval sin 2 #i3~ 10~ 5 — 10~ 2 [38, 39]. Here 612 — 34° is the 
"solar" mixing angle (see e.g. [42]). For the extreme values of #13, the dependence of 
the survival probability on the energy is negligible in first approximation; oscillations in 
the Earth are also at the level of few per cent when averaged over the different arrival 
directions of the neutrinos [43], therefore they will be neglected here. 

For Am| x < the calculation of p is complicated by the effects of neutrino -neutrino 
scattering, which could be strong few seconds after the start of the neutrino burst [44] 
and induce a swap of the v e and v x fluxes above a critical energy E c ~ 6 — 9 MeV 
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[30, 31, 32, 33, 34, 35, 36]. Still, the survival probability remains within the interval in 
Eq. (3.3) at all times (see e.g. [34]), so one can always define a time averaged survival 
probability with value in this interval. 

For our purposes, it is sufficient to consider p as energy independent, as E c is 
typically below threshold of the cross sections of interest (see sec. 4), and take the two 
extreme cases of total (p = 0) and minimal (p = sin 2 # 12 ) flux permutation. This is 
adequate to show the extent of variation of the 97 Tc abundance with the conversion 
pattern. 



10 b 



10 4 



T 
> 

d) 100: 



I 



E 
o 



0.01 




Figure 2: Time averaged solar (leftmost curve) and supernova (lines extending to the right) 
u e fluxes at Earth, inclusive of oscillation effects. The solid curves correspond, from upper to 
lower, to the best, natural and worst scenarios in Table 2 and complete flavor permutation 
{p = 0, see Eqs. (3.2) and (3.3)). The dashed ones refer to the same scenarios with the 
minimal permutation (p ~ 0.31). A rate R sn = 10~ 2 yr" 1 was used. 

Fig. 2 shows the time averaged supernova v e flux at Earth: 

Fr th = R sn F e /(A7Td 2 eff ) . (3.4) 
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We have plotted the six fluxes at Earth obtained with the two extreme values of p 
and the three scenarios in Table 2. They are averaged over time for the purpose of 
comparison with the solar neutrino flux 1 . 

All that was discussed so far refers to the most common scenario of core collapse 
supernova, the one that leads to a successful explosion and the formation of a neutron 
star. Stars of mass above 25-40 M (M = 1.99 • 10 30 Kg being the mass of the sun), a 
10 — 20% fraction of all collapse candidates, could either explode and form a black hole 
by fallback or collapse into a black hole directly, with no explosion (see e.g. [45] and 
references therein). The neutrino emission in the latter case has been studied recently 
[46, 47, 48, 49, 50] and found to be characterized by a much shorter burst (0(1) s or 
less) with luminosity and average energies higher than the neutron-star forming case, 
especially in v e and v e . Higher and more energetic fluxes have been found [46, 47, 48, 49] 
for the stiffer equation of state (EoS) by Shen et al. (S) [51] compared to the softer 
one by Lattimer and Swesty (LS) [52]. 

Here we model the contribution of these jailed supernovae to the time averaged 
supernova v e flux (shown in fig. 3) using the original fluxes presented in ref. [50] 
(fig. 5 there). Those have the character of examples, only roughly representing the 
possible variations that one can have with the EoS and progenitor model, for which a 
comprehensive study still lacks. They are characterized by the parameters Eq x ~ 24 
MeV, L x ~ 0.45 • 10 53 ergs, E 0e ~ 20.7 MeV, L e ~ 1.4 • 10 53 ergs for the S EoS and by 
E 0x ~ 22 MeV, L x ~ 0.2 ■ 10 53 ergs, E 0e ~ 17 MeV, L e ~ 0.54 • 10 53 ergs for the LS 
one. For the spectral shapes, we checked that Eq. (3.1) is a reasonable approximation; 
still we chose to interpolate the numerically calculated points in [50] with a polynome 
of order 4 in the plane log(-E') — log(F°). Results have only minor differences compared 
to the more conservative linear fit used in [50]. 

In absence of indications of the opposite, we assumed that black hole forming 
collapses have the same spatial distribution, Eq. (2.1), as the neutron star forming ones, 
and that the estimated ratio of rates for the two types of collapses, Rbh/R sn ~ 0.1 — 0.2 
is valid for our galaxy. We have also considered the conversion effects to be the same 
for the two types, Eqs. (3.2) and (3.3). While this is not the case for every value of 
#13, it is true for the extreme values of #13 and of p, which are of the most interest 
here [50]. Fig. 3 reveals that the contribution of failed supernovae is at the level of 
tens of per cent of the total flux at E ~ 20 — 30 MeV, and can exceed the flux from 
regular supernovae above 50-60 MeV depending on the parameters. Notice that for 

^^Using a time averaged flux might be somewhat confusing, because it might suggest the idea of 
a flux that is continous in time. In constrast with the solar one, the supernova neutrino flux we are 
considering is far from continous: it consists of ~ 10 s bursts that reach the Earth a few times per 
century and whose effects on the 97 Tc production accumulate over several millions of years. 
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Figure 3: The same of fig. 2 for the flux of failed supernovae, obtained using the fluxes in 
fig. 5 of [50]. Of the two solid curves, the upper (lower) is for the Shen (Lattimer-Swesty) EoS. 
The same holds for the dashed lines. We took a rate for failed collapse of Rbh = 10 -3 yr . 

failed supernovae the v e flux is smaller for total flavor permutation p = 0, because in 
this case the flux of u e s at Earth originates from the less luminous original v x flux. 

3.2 Solar neutrinos 

Solar neutrinos are the only significant background that can not be eliminated. In 
the energy range of interest here, z/ e s from the 8 B and hep reactions of the fusion 
chain are relevant. The former component extends to 16.6 MeV energy, while hep 
neutrinos reach 18.8 MeV. We describe these fluxes following the Bahcall-Serenelli- 
Basu BS05(PO) solar model [53]. There, the 8 B component is normalized to = 
5.69 • 10 6 cm~ 2 s _1 , and the hep one to $hep — 7.93 ■ 10 3 cm~ 2 s _1 (unoscillated values). 
The solar neutrino survival probability as a function of energy was taken as the best fit 
in [54] (fig. 2 of this reference), corresponding to the LMA-MSW oscillation parameters 
Am^ = 7.92- lCr 5 eV 2 and sin 2 12 = 0.314. With these, the oscillated fluxes at Earth 
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are $^ sc = 2.17 • 10 6 cm^s" 1 , and $^ = 2.81 • 10 3 cm^s" 1 . The total of the two, 
differential in energy, is shown in figs. 2 and 3. 

Notice that the oscillated fluxes at Earth used here differ substantially from those 
used by Haxton & Johnson: $f J = 1.1 • 10 6 cm~ 2 s _1 , and ^ e J p = 1.6 • 10 4 cm~ 2 s~ 1 . 
This difference is a result of the two decades of progress in solar neutrino research, and 
is one of the major novelties of our work. 

A comment on uncertainties on the solar neutrino flux is due. Errors on the mea- 
sured Amj! and sin 2 6*12 are of ~ 3% and ~ 10% (at la) respectively [42]; corresponding 
to an error of the order of ~ 10% on the solar neutrino survival probability. The solar 
model itself has uncertainties, e.g. those associated to using different inputs for the 
solar opacities [53], which are at the level of (9(0.1)%. Here we do not include all 
these uncertainties explicitly, because they do not impact our conclusions significantly; 
moreover, we note that for our purposes the oscillated fluxes $^ and are the 
only relevant quantities (assuming their spectra as known). These are measurable very 
precisely, in principle. The SNO experiment has already measured (the part of the 
flux above its threshold of ~ 5 MeV) with precision of ~ 3-5% (statistical and system- 
atic, la) [42], while a measurement of hep neutrinos requires to wait for the next phase 
of solar detectors [55]. The importance of a precise measurement of the solar neutrino 
fluxes in view of a possible 97 Tc experiment will be discussed briefly in sec. 5. 

4. Neutrino-induced cross sections 

Predictions of neutrino-nucleus cross sections in the solar and supernova energy range 
require a precise knowledge of the transitions to the low lying states and to the res- 
onance region, well above the particle emission thresholds. The latter include both 
allowed transitions associated to the Gamow- Teller and Isobaric Analogue State - and 
the forbidden ones, such as the spin-dipole and higher multipoles. The particular reac- 
tion cross sections of interest are 98 Mo(z/ e , e~n) 97 Tc and 97 Mo(u e , e~) 97 Tc. In this work 
we follow two different procedures to calculate the contribution to the cross sections 
coming from the Gamow- Teller and Isobaric Analogue State on one hand and the spin- 
dipole and higher multipoles (up to J = 5) on the other hand. In the first one we use 
the allowed approximation, which corresponds to the low momentum transfer, as done 
in Haxton and Johnson's paper. For all other multipoles we use the Walecka formalism 
[56] for the reaction cross section without neglecting the momentum transfer, with the 
transition matrix elements provided by the microscopic Quasi-Particle Random-Phase 
Approximation (QRPA). Note that this is at variance with [1] where only the spin- 
dipole is included whose contribution is estimated using the Goldhaber- Teller model. 
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Figure 4: 98 Mo(f e , e~n) 97 Tc reaction cross section as a function of the neutrino energy. 
The dashed curve presents the contribution due to the Gamow- Teller and Isobaric Analogue 
State transitions, evaluated within the allowed approximation using experimental data. The 
dotted curve shows the contribution coming from the spin-dipole and other multipoles (up 
to J = 5), calculated with the QRPA approach. The total cross section is also shown (solid 
line) . 

To evaluate the transition matrix elements of the Gamow- Teller and Isobaric Ana- 
logue State within the allowed approximation, we exploit the only measurement avail- 
able, the one by J. Rapaport et al. [57], who studied 98 Mo(p, n) 98 Tc reaction at forward 
angles for 120 MeV incident protons. This procedure is the same as in [1]. It is im- 
portant to remind that the model-independent Fermi and Ikeda sum-rules tell us that 
the total strength associated to the Fermi and Gamow- Teller transitions is given by 
N-Z and 3(N-Z) respectively, N and Z being the number of neutrons and protons of 
the considered nucleus. In [57] it is found that a fraction of the Gamow- Teller strength 
B(GT), i.e. ^B(GT)=28 ± 5, is distributed to excite the 98 Tc states having up to 18 
MeV above the ground state. Since the experimental data are cut above 18 MeV we 
have checked the influence of a B(GT) strength located at higher energy by extrapolat- 
ing the experimental data with different curves (Gaussian or Lorentzian). The results 
(sec. 5) turn out not to be sensitive to the specific extrapolation curve chosen. On the 
other hand, the Fermi transitions are dominated by the Isobaric Analog State of 98 Mo, 
located at E = 9.656 MeV above 98 Tc ground state. This single resonance carries all 
the Fermi strength B(F) given by the Fermi sum-rule, namely B(F) = (N-Z) = 14. 
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Concerning the forbidden transitions, their contribution to the cross section has been 
calculated using the QRPA approach. The largest contribution comes from the 1~ and 
T~ transitions to 98 Tc which essentially affects the high energy part of the cross section 
(above 40 MeV). The total cross section we use in this work is shown in Figure 4. 

The need to minimize the background from solar neutrinos motivates the choice 
of the process 98 Mo(z/ e , e~n) 97 Tc, where one neutron is emitted. For this reaction the 
threshold in energy for neutron emission is high enough (8.96 MeV) to suppress the 
pp solar neutrino flux completely and the 8 B and hep fluxes partially compared to the 
supernova contribution. We make the assumption that all the excited states of 98 Tc 
laying above neutron emission threshold will decay by emitting neutrons and produce 
97 Tc in its ground state. Note that the 97 Mo+p threshold is by 1.1 MeV lower than 
97 Tc+n one. However proton emission at low energies should be highly suppressed 
compared to the neutron one due to the Coulomb barrier. At higher energies the 
97 Mo+p channel should become as much important as the 97 Tc+n one. Thus, at these 
energies, our predictions should overestimate the rate of 97 Tc production from 98 Mo. 

A delicate task is to accurately predict the 97 Mo(z/ e , e~) 97 Tc cross section, since no 
experimental information is available. Following ref [1] for this case we have evaluated 
the 97 Mo(i/ e ,e") 97 Tc cross section exploiting the sum-rule argument, namely scaling 
the B(GT) strength from the 98 Mo(p, n) 98 Tc measurement by the ratio (N-Z)97 Mo /(N- 
Z)98 Mo . This scaling is performed after having subtracted the near-threshold events 
describing the transition to to the 98 Tc ground state. In fact, the 97 Mo(z/ e , e~) 97 Tc 9 . s . 
transition matrix elements is negligibly small as one can check by converting the 97 Tc 9 . s . 
lifetime {t\/2 = 2.6 • 10 6 s). Concerning the Fermi strength, we consider that all the 
strength B(F)=2T=13 is associated to the Isobaric Analogue State located at 11.05 
MeV excitation energy of the 97 Tc nucleus. As in the 98 Mo case, we neglect 97 Tc 
production above the 96 Tc+n threshold. Since again we make the assumption that 
100% of the cases a neutron is produced, this time we slightly underestimate the 97 Tc 
production from 97 Mo, which will give in particular fewer supernovae events. 



5. Results 

5.1 Rates 

The table in figure 5 presents our predicted 97 Tc production rates from 98 Mo and 97 Mo 
due to both solar and supernova neutrinos. 

Let us first comment on the 98 Mo(z/ e , e~n) 97 Tc channel. Solar neutrinos produce 
97 Tc essentially via the Gamow- Teller and the Fermi transitions to excited states of 
98 Tc, which then decay by emitting a neutron; the contribution from forbidden tran- 



- 12 - 





98 Mo(zv e ,e-n) 97 Tc 


97 Mo(i/ e ,e-) 97 Tc 




HT 


Allowed 
only 


Total 

(allowed + 
forbidden) 


HT 


Total 

(allowed + 
forbidden) 


8 B 


3.17 


5.07 


5.07 




37.4 


hep 


1.0 


0.14 


0.14 




0.2 


Total solar 


4.17 


5.21 


5.21 


31 


37.6 


SN (best) 




0.73 - 1.0 


1.13 - 1.58 




0.23 - 0.30 


SN (natural) 


5.4 


0.28-0.36 


0.37 - 0.49 


1.7 


0.10-0.12 


SN (worst) 




0.12-0.14 


0.15-0.18 




0.048 -0.051 


BH (S Eos) 




0.08 -0.12 


0.17-0.22 




0.02 - 0.03 


BH (LS EoS) 




0.03 - 0.04 


0.06 - 0.07 




0.08-0.010 



Figure 5: 97 Tc production rates (10 -37 x s" 1 ) associated to neutrinos from the sun and 
supernovae (where BH indicates failed, black hole- forming supernovae), in comparison with 
Haxton and Johnson's (HJ) ones. Intervals for SN rates take into account the variations in the 
neutrino oscillation pattern (normal vs inverted mass hierarchy and small vs large #13). All 
our results are obtained with a galactic rate R sn = 0.01 yr" 1 (Rbh = O.OOlyr" 1 ) for successful 
(failed) supernovae. They can be easily rescaled for other values of this quantity. Natural 
values of R sn are 0.02 - 0.03 yr" 1 , but the larger interval 0.01 - 0.05 yr 1 is allowed (table 
1). HJ's numbers refer to a supernova rate of 0.09 yr -1 . 

sitions is negligible. In comparison to HJ [1], we find a higher rate of transition due 
to 8 B neutrinos, result of the larger (almost double) normalization of this flux after 
oscillations (sec. 3). In contrary, the rate due to the hep flux is an order of magnitude 
smaller because of the smaller initial normalization and the further reduction of the 
flux due to flavor conversion. 

In Figure 5 we show predictions for supernova neutrinos considering a galactic 
supernova rate of 0.01 yr" 1 and taking into account a range of values for the supernova 
neutrino fluxes at the neutrino-sphere and at the surface of the supernova, due to the 



13 



unknown neutrino mixing parameters (Table 2, fig. 2). Contrary to the solar case, the 
rate from supernova neutrinos has an important contribution (~30-60%) coming from 
the forbidden transitions since the neutrino spectra cover several tens of MeV energy 
range. Obviously, an increased luminosity or average neutrino energy boosts the 97 Tc 
production rate. One can see that our predictions are typically one order of magnitude 
smaller than what Haxton and Johnson estimated. The difference is mostly due to a 
higher supernova rate, R sn = 0.09 yr _1 and a higher neutrino luminosity (8 • 10 52 ergs 
compared to our Lq = 0.5 ■ 10 53 ergs), as well as the fact that the flavor permutation 
of supernova neutrinos was not included by HJ. Inclusion of oscillations increases the 
supernova flux above ~ 15 — 20 MeV, as observed in sec. 3, and therefore lessens the 
gap between our results and those in [1] 2 . Still, while HJ concluded that the rate of 
98 Mo(^ e , e~n) 97 Tc reaction due to supernova neutrinos should be as much as 1.3 times 
the solar one, our finding is that the supernova contribution is most likely at the level 
of tens of per cent, but can amount to more than half the solar contribution accepting 
optimistic neutrino parameters and R sn ^ 0.03 yr _1 , as allowed by current estimates 
(table 1). 

The effect of failed supernovae on the 9S Mo(u e , e~n) 97 Tc rate is significant. While 
certainly small compared to the solar rate, for the S EoS this contribution can easily 
reach ~ 50% of the one from successful supernovae, and can be even of the same size 
for the largest fraction of failed versus successful supernovae, Rbh/ l Rsn ~ 

0.2. For the 

LS EoS the effect is more modest, but still at the level of ~ 10 — 30% of the total 
supernova contribution. We notice that, due to the more energetic neutrino fluxes, the 
rate from failed supernovae is very sensitive to the forbidden states, and nearly doubles 
when these are included in the cross section calculation. 

The 97 Tc production via the 97 Mo(z/ e , e^) 97 Tc channel constitutes a very different 
case, also illustrated in fig. 5. For this process the reaction threshold is only 320 keV, 
and therefore a large fraction of the solar 8 B and hep neutrino fluxes contribute to it. 
Moreover at 9.8 MeV the 97 Mo(i/ e , e"n) 96 Tc channel opens and starts dominating over 
the 97 Tc production, further suppressing the supernova contribution to this channel. 
For these reasons the 97 Tc production from 97 Mo is largely dominated by solar neutri- 
nos, primarily due to the GT transitions, with the supernova-induced production being 
only a few per cent of the total. The effect of the 97 Mo(z/ e , e~) 97 Tc channel is mostly 
to dilute the supernova signal due to scattering on 98 Mo into a larger background and 
make a measurement more challenging. 



2 If we make the same assumptions as in [1] we get the same solar background (with differences 
up to 8%) and a supernova rate about 17% lower. This is probably due to the different procedure 
followed in the present work as far as the forbidden transitions are concerned. 
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5.2 Abundances of 97 Tc in Molybdenum rocks 

In order to estimate the feasibility of a geochemical experiment one should translate 
the calculated rates into the number of 97 Tc atoms in a sample of rock, which is the 
experimentally measurable quantity. We calculated this for 10 kt mass of Molybdenum 
ore, considering that it contains about 49 t of M0S2 [1], and assuming that its Mo 
content follows the natural isotopic abundances (i.e., 24.13% for 98 Mo and 9.55% for 
97 Mo). 





98 Mo(i/ e ,e-n) 97 Tc 


97 Mo(*/ e ,e-) 97 Tc 


Total 


HJ, 
total 


Solar only 


2.7 


7.8 


10.6 


8.0* 


Solar 
+ 100% SN 


3.5 


7.9 


11.4 


11.2 


Solar + 90% SN 
+ 10% BH (S 
EoS) 


3.7 


7.9 


11.6 





* 8 B only. See [1] for details. 



Figure 6: Numbers of 97 Tc atoms (10 6 unit) expected in 10 kt of rock from the two 
different processes of interest. Totals are given too, as well as the results for the capture 
rates of Haxton and Johnson (HJ) for comparison. Our results refer to a supernova rate of 
0.03 yr" 1 and parameters in the "natural" scenario (Table 2) for successful supernovae (SN). 
Results in the last row refer to a scenario with a 10% failed supernovae (BH) with neutrino 
emission as in the S model of ref. [50] (see fig. 3). In all cases we have considered complete 
flavor permutation (p = 0), which gives the largest 97 Tc production rates for successful SNe. 

Results are presented in fig. 6. It appears that in 10 kt of Molybdenum ore that 
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is completely shielded from cosmic rays, there should be ~ 1.1 • 10 atoms of Tc, 
of which about 10% from supernovae, for natural parameters of supernova neutrino 
emission and R sn = 3 • 10~ 2 yr _1 . This contribution can be larger by up to a factor 
of 3 if the best scenario of supernova neutrino spectra and luminosities (table 2) is 
realized. The effect of failed supernovae is a few per cent of the total abundance. We 
note how our calculated abundance is very close to that of Haxton and Johnson: this 
similarity is completely accidental, since our input quantities are very different and we 
find a substantially higher solar contribution, balanced by a lower supernova one. 

Let us now discuss the significance of our result in view of an experiment. The 
first question is if the neutrino-induced 97 Tc is at detectable level. The answer is 
affirmative: already at the time of HJ the best instruments were quoted [1] to have 
sensitivity down to about 5 • 10 6 atoms for a pure Tc sample. A modern experiment 
would employ essentially the same methods, with a sensitivity slightly improved (a 
factor of two or so) by a number of technological updates [58]. Therefore we come 
to the same conclusion as HJ: at the price of setting up an industrial-scale extraction 
project, capable to process about 10 kt of rock at large depth 3 , an experiment should 
be feasible. Such experiment should be arranged in such a way to avoid the "roast 
memory" problem that was fatal to the original experimental attempt in 1988 [5]. 
While potentially laborious, this is not difficult conceptually. 

Once established that a detection is possible, a key question is if the supernova 
neutrino contribution can be unambiguously identified. Our study makes it clear that 
this contribution is comparable to the several uncertainties involved in the problem. In- 
deed, already the GT strength distribution obtained from 98 Mo(p, n) 98 Tc measurement 
has more than 17% uncertainty, whereas the 97 Mo(V, e~) 97 Tc cross section estimate 
along with branching ratios into the decay channels is purely phenomeno logical. One 
should also consider the ~ 5 — 10% error on the oscillated 8 B solar neutrino flux, and 
the still unmeasured hep flux. 

Therefore, we conclude that a measurement of the 97 Tc abundance in deep Molyb- 
denum ore would be an interesting complement to studies of the solar neutrino flux, 
but would provide only loose upper limits on the supernova flux, perhaps at the level 
of the constraint from neutrino detectors on the supernova rate (table 1). Two major 
experimental steps should be achieved in the near future to render the geochemical 
quest for galactic supernova neutrinos attractive : i) an increase in precision in the 
measured 8 B solar neutrino flux and a new precise measurement of the hep flux; ii) a 
precise meaurement of the transition matrix elements for the Gamow- Teller as well as 
the forbidden multipoles contributing to the 97 Mo(p, n) 97 Tc and 98 Mo(p, n) 98 Tc reac- 

3 We stress that the use of industrial-scale equipment and large scale excavations - while exceeding 
the typical scale of geophysics projects - are common in the field of neutrino astrophysics. 
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tions, along with the branching ratio relevant to 97 Tc production. This can be obtained 
from future charge-exchange measurements for example at KVI or RIKEN. A direct 
measurement of the neutrino induced cross sections could also be performed at future 
low energy neutrino facilities exploiting conventional sources, such as at z/SNS, or at a 
low energy beta-beam facility [59, 60, 61]. 

6. Conclusions 

In this paper we have investigated the possibility to observe the integrated supernova 
galactic neutrino contribution by measuring the 97 Tc produced by electron neutrinos in 
a Molybdenum ore, as first proposed in [1]. With this aim realistic estimates - reflecting 
the state of the art in the neutrino physics and of galactic supernovae neutrino fluxes - 
have been performed. For our calculations we have used experimental information on 
the Gamow- Teller and Fermi distributions of the 98 Mo(z/ e , e~) 98 Tc reaction, combined 
with QRPA predictions for the forbidden contributions. 

We have found that the relative production of 97 Tc from supernova neutrinos is 
about 0.1 of the solar background, which is much smaller than what was evaluated 
in the pioneering study [1]. For a geochemical measurement to become attractive, 
a better knowledge of the solar neutrino spectrum and of the 97 Mo(^ e , e~) 97 Tc and 
9S Mo(u e , e~n) 97 Tc cross sections are necessary. An improved experimental sensitivity 
is desirable to reduce the quantity of rock to be processed. Future precise charge- 
exchange measurements at high energies and forward angles and/or neutrino-nucleus 
measurements for the 97 Mo and 98 Mo isotopes will defintely be very helpful. 

We are indebted to W. Haxton and A. Anbar for very fruitful exchanges. C.V. and 
R.L. acknowledge the support from "Non standard neutrino properties and their impact 
in astrophysics and cosmology", Project No. ANR-05-JCJC-0023. CL. acknowledges 
ASU and RBRC for support, and is grateful to the IPN (Orsay) for generous hospitality 
while working on this project. 
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